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An effective field theory for the three-body system with large scattering length 
is applied to three-body recombination to a weakly-bound s-wave state in a 
Bose gas. Our model independent analysis demonstrates that the three-body 



■ recombination constant a is not universal, but can take any value between 

zero and 67.9 ha 4 /m, where a is the scattering length. Other low-energy three- 
{Sj , body observables can be predicted in terms of a and a. Near a Feshbach 

resonance, a should oscillate between those limits as the magnetic field B 

o 

approaches the point where a — > oo. In any interval of B over which a 
' increases by a factor of 22.7, a should have a zero. 

a ■ 
g. 

The successful achievement of Bose-Einstein condensation has triggered a large interest 
in interacting Bose gases. One of the main factors limiting the achievable density in these 
experiments is the loss of atoms through 3-body recombination. Such losses occur when 
three atoms scatter to form a molecular bound state (called a "dimer" for brevity) and a 
third atom. The kinetic energy of the final state particles allows them to escape from the 
trapping potential. This 3-body recombination process is interesting in its own right as it 
provides a unique window on 3-body dynamics. 

The number of recombination events per unit time and volume can be parametrized 
as v rec = an 3 , where a is the recombination constant and n the density of the gas. The 
calculation of a in general is a complicated problem, because it is sensitive to the detailed 
behavior of the interaction potential The simplest case is 3-body recombination to a 
weakly-bound s-wave state. For atoms of scattering length a and mass m, the binding energy 
of the dimer is = h 2 / (ma 2 ) and the size of the bound state is comparable to a. Assuming 
that a is the only important length scale, dimensional analysis implies a = Cha 4 /m, where 
C is dimensionless. The problem of 3-body recombination to a weakly-bound s-wave state 
has been studied previously Fedichev et al. found that the coefficient C has the 

universal value C = 3.9. Nielsen and Macek [BJ and Esry et al. M found that C could take 
any value between and C max ^65. 

An ideal means to study the dependence of a on the scattering length a is to use Feshbach 
resonances |J, which occur when the energy of a spin-singlet molecular bound state is tuned 
to the energy threshold for two separated atoms by applying an external magnetic field B. 
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Such resonances have, e.g., been observed for 23 Na and 85 Rb atoms 0f|. When the magnetic 
field is varied in the vicinity of the resonance, the scattering length varies according to || 

a(B) = ao (l + , (1) 

where a is the off-resonant scattering length and A and B characterize the width and 
position of the resonance, respectively. On the side of the resonance where a increases 
towards +00, the spin-singlet molecule becomes a weakly-bound s-wave state. 

In this Letter, we use effective field theory methods to calculate the recombination con- 
stant a for a weakly-bound s-wave state. We find that the naive scaling relation a = Cha 4 / m 
is modified by renormalization effects, so the coefficient C is not universal, and that its maxi- 
mum value is C max = 67.9. Our results confirm those of Refs. |||4j], demonstrate that they are 
model-independent, and provide a new tool for precise studies of the recombination rate. On 
the side of a Feshbach resonance where a — > 00, the -B-dependence of a can be predicted in 
terms of one adjustable parameter. It has a remarkable behavior, with C oscillating between 
zero and 67.9 more and more rapidly as B approaches the resonance. 

Effective field theory (EFT) is a powerful method for describing systems composed of 
particles with wave number k much smaller than the inverse of the characteristic range R 
of their interaction. (For a van der Waals potential -C 6 /r 6 , R = (2mC 6 /ft 2 ) 1 / 4 ). EFT 
focuses on the aspects of the problem that are universal, independent of the details of 
short- distance interactions, by modelling the interactions as pointlike. The separation of 
scales < allows a systematic expansion in powers of the small parameter kR ||. 
Generically, the scattering length a is comparable to R, and the expansion is effectively in 
powers of ka. The pointlike interactions of the EFT generate ultraviolet divergences, but 
they can be absorbed into the renormalized coupling constants of the effective Lagrangian. 
All information about the influence of short- distance physics on low-energy observables is 
captured by these constants. At any given order in kR, only a finite number of coupling 
constants enter and this gives the EFT its predictive power. The domain of validity of EFT is 
kR <C 1, even in the case of large scattering length a> R. Thus it should accurately describe 
weakly-bound states with size of order a. However, the dependence on ka is nonperturbative 
for k ~ 1 1 a, and it is necessary to reorganize the perturbative expansion into a new expansion 
in kR by resumming higher order terms to all orders in ka. There has been significant 
progress recently in carrying out this resummation for the 3-body system. At leading order 
in kR, a single 3-body parameter is necessary and sufficient to carry out the renormalization 
Til . The scattering length and this 3-body parameter are sufficient to describe all low- 



energy 3-body observables up to errors of order R/a. In nuclear physics, this result has 
allowed a successful description of low-energy neutron-deuteron scattering and the binding 
energy of the triton. The variation of the 3-body parameter provides a natural explanation 



for the Phillips line [|TT 



We will apply this EFT for systems with large scattering length to the 3-body recombi- 
nation problem. For simplicity, we now set h — 1. We start by writing down a general local 
Lagrangian for a non-relativistic boson field ip: 
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FIG. 1. First row: diagrams contributing to 3-body recombination. Second row: integral 
equation summing them to all orders. Atom propagator and exact dimer propagator are indicated 
by thin (thick) lines, respectively. 

The dots denote terms with more derivatives and/or fields; those with more fields will not 
contribute to the 3-body amplitude, while those with more derivatives are suppressed at low 
momentum. In order to set up integral equations for 3-body amplitudes, it is convenient to 
rewrite L by introducing a dummy field d with the quantum numbers of two bosons, 

C = ^(i^- + X^U + Sd - -4=(dtyV + h-c.) 
V at 2mJ v2 

+ h<#diflij) + ... . (3) 

The original Lagrangian (fj) is easily recovered by a Gaussian path integration over the d 
field, which implies d = (g / \fi)i\) 2 / (\ + hip^ip), C = g 2 , and D = —3hg 2 . The atom 
propagator has the usual non-relativistic form i/(uo — p 2 /2m). The bare propagator for d is 



simply i, but the exact propagator, including atom loops to all orders, is |10| 



■ Q , -iAit/(mg 2 ) 

tS d {uj,p) = = = , (4) 

— 1/a + J—muj +p 2 /A — ie 

where a is the scattering length, which is related to the bare parameter g and the ultraviolet 
(UV) cutoff A by 



mg / mg A 
ItT V + 2tt 2 



(5) 



The propagator @ has a pole at u — —1/ (ma 2 ) +p 2 / (Am) corresponding to a weakly-bound 
state. Attaching four atom lines to this propagator gives the exact two-atom scattering 
amplitude. Taking the incoming atoms to have momenta ±p, the amplitude is (— 1/a— ip)^ 1 , 
confirming the identification of a as the scattering length. 

We now consider the 3-body recombination process. We take the momenta of the incom- 
ing atoms to be small compared to the momenta of the final particles, which have magnitude 
Pf. Using Fermi's golden rule, the recombination coefficient can be written 

a = (6) 
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where T(p) is the amplitude for the transition between three atoms at rest and a final state 
consisting of an atom and a dimer in an s-wave state with momentum p in their center-of- 
momentum frame. In Eq. (|6]), this amplitude is evaluated on shell at the value Pf = 2/(y / 3o) 
prescribed by energy conservation. However, T(p) is also defined at off-shell values of p. The 
first few diagrams contributing to T are illustrated in Fig. [I]. All loop diagrams are of the 
same order as the tree diagrams, and they therefore have to be summed to all orders. This 
is conveniently accomplished by solving the integral equation represented by the second 
equality in Fig. |l|. Corrections to this equation are of order R/a. The integral equation is 
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q 2 + pq + p 2 



q 2 — qp + p 2 



+ 



h 



where we have inserted an UV cutoff A on the integral over q. If we were allowed to 
set h = and take A — > oo, a rescaling of the variables in Eq. fl7|) would lead to T(jp) = 
K (ap)a 3 / 2 1 (mp) , with K{x) a dimensionless function. Evaluating this in Eq. ([]), the scaling 
relation a = Ca 4 /m would follow immediately. However, the integral equation (^) has the 
same properties as the one describing atom-dimer scattering [10] and the limit A — > oo can 
not be taken. The individual diagrams are finite as A — > oo, but their sum is sensitive to 
the cutoff. In an EFT, the dependence on the UV cutoff is cancelled by local counterterms. 
In Ref. [|TIJ, it was shown that the dependence of the low-energy observables on the cutoff 
A could be precisely compensated by varying h appropriately. Writing h = 2mg 2 H(A)/A 2 , 



it was found that H(A) could be well approximated by 

H(A) w - tan [s ln(A/A*) - tt/4] 



(8) 



where sq ~ 1.0064 is determined by the asymptotic behavior of the integral equation. This 
expression defines a parameter A* |K| that characterizes the effect of the 3-body force on 
physical observables. A remarkable feature of this expression is its periodicity in In A. As A is 
increased, H(A) decreases to — oo, changes discontinously to +oo, and continues decreasing. 

The scaling violations in T(p) from the renormalization of the 3-body force can be ex- 
pressed as a dependence on aA*. The simple scaling relation for a is therefore replaced by 
a = C(A*a) a 4 /m. Thus the value of C is not universal. In Fig. ^|, we show a as a function 
of a for A*ao = 1.78, 4.15, 7.26 and 19.77, where is an arbitrary but fixed length scale. 
Interestingly, a appears to oscillate as a function of In a between zero and a maximum value 
C max - We find that the curves can be very well fit by the expression 



a 



ha 



rn 



C max cos 2 [s ln(aA*) + 5] , 



(9) 



with C max = 67.9 ± 0.7 and 5 = 1.74 ± 0.02. 

We now compare our result to those of Refs. Using an approximate solution to 

the 3-body wavefunction in hyperspherical coordinates, Fedichev et al. found that the 
coefficient C has the universal value C = 3.9, independent of the interaction potential. We 
find that the value of C is not universal, but can vary from zero to about 67.9, depending 
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FIG. 2. a as a function of a for different values of the 3-body parameter: clqA* = 1.78 (full 
line), 4.15 (dotted line), 7.26 (dashed line), and 19.77 (long-dashed line). 

on the details of 3-body interactions at short distances. We suggest that the specific value 
in Ref. @ must correspond to an implicit assumption about the short-distance behavior of 
the 3-body wavefunction. Our results are consistent with those of Refs. [@,£|, which were 
obtained using the hyperspherical adiabatic approximation. Nielsen and Macek |J obtained 
C max ~ 68.4 by applying the hidden crossing theory. Esry et al. J|] used coupled channel 
calculations to obtain a numerically for over 120 different 2-body potentials. For a ^> R, 
their empirical result has the form of Eq. @ with sq — 1 and C max = 60 ± 13. Refs. 

show that the zeroes in a arise from interference effects involving 2-body and 3-body 
hyperspherical adiabatic potentials. The origin of these interference effects is less obvious 
in our EFT approach. However, our approach has several other advantages. First, it is 
completely model independent. Second, it is a controlled approximation, with corrections 
from finite range effects suppressed by powers of R/a. Third, it has predictive power in 
that all other low-energy 3-body observables can be determined in terms of a and A*. For 
example, the atom-dimer scattering length a d can be fit by JTOj 

a d ^a (1.4 - 1.8 tan[s ln(aA*) + 3.2]) . (10) 

We now apply the EFT to Feshbach resonances, where the value of a is varied by changing 
the external magnetic field (cf. Eq. ([I])). Our formalism is only applicable close to the 
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Feshbach resonance on the side where a > 0, because only in that region is there a weakly- 
bound molecule with ~ l/(ma 2 ). Away from the resonance, or close to the resonance but 
on the side where a < 0, 3-body recombination must involve more deeply-bound molecules 
with binding energies of order l/(mR 2 ). To predict a as a function of the magnetic field 
B, we must specify how the parameter A* in Eq. (|8|) varies as a function of B. A Feshbach 
resonance is characterized by nonanalytic dependence of the scattering length a on B. In 
a field theory, nonanalytic dependence on external parameters arises from long-distance 
fluctuations [12|]. The explicit ultraviolet cutoff A in our EFT excludes long-distance effects, 
which could introduce nonanalytic dependence on B, from the coefficients in C Thus the 
bare parameters Co = g 2 and D = —3hg 2 in Eq. (0) should be smooth functions of B for a 
fixed value of A. The resonant behavior of the scattering length near B = Bq in Eq. (|l]) can 
be reproduced by approximating g 2 by a linear function of B in the resonance region. The 
parameter h should also be a smooth function of B, but we can take h to be approximately 
constant over the narrow resonance region. This assumption implies via Eq. (||) that A* 
should be constant while a(B) varies as in Eq. (||). 

The behavior of the recombination coefficient near the Feshbach resonance can therefore 
be read off from Fig. 0. If a is measured at one value of B for which a(B) ^> |ao|, it 
determines a as a function of B up to a two-fold ambiguity corresponding to whether the 
slope of C is positive or negative at that value of B. As B approaches B , C should oscillate 
between and C max in the manner shown in Fig. [2|. The successive zeros correspond to values 
of a(B) that differ roughly by multiplicative factors of exp(7r/so) ~ 22.7. Thus EFT makes 
the remarkable prediction that there are values of the magnetic field close to a Feshbach 
resonance where the contribution to a from a weakly-bound state vanishes. 

The loss rate of 23 Na atoms from a Bose-Einstein condensate near a Feshbach resonance 
has been studied by Stenger et al. || . Our theory applies to the low- field side of the resonance 
at B = 907 G. Taking into account the 3 atoms lost per recombination event and a Bose 
factor of 1/3!, the loss rate from the condensate is N = —aN(n 2 )/2. The loss rates measured 
in Ref. || correspond to a coefficient C ~ 300 both off and near the resonance. This value 
is a factor of 4 larger than our maximum value. If C > C max , the 3-body recombination rate 
must be dominated not by the weakly-bound Feshbach resonance but instead by molecules 
with much larger binding energies ~ l/(mR 2 ). Alternatively, the large loss rate in Ref. || 
could be due to collective effects associated with the Bose-Einstein condensate, such as a 

2- body recombination process involving atoms and dimers in a molecular condensate ||. In 
Fig. |3], we show the contribution to a from the weakly-bound state as a function of the 
magnetic field B for A* = 1.94 nm" 1 , 1.16 nm -1 , and 0.65 nm" 1 . If this contribution to a 
could be isolated, the first zero may be wide enough to be observed by experiment. The 
higher zeroes, however, are increasingly narrow and very close to the resonance. 

We have applied an EFT for atoms with large scattering length to the problem of 3- 
body recombination into a weakly-bound s-wave state. We find that the coefficient C in the 
scaling relation a = Cha 4 /m is not universal, but must be in the range < C < 67.9. If the 

3- body recombination rate is measured to be larger than the maximum value, there must be 
a large contribution from molecules that are more deeply bound. Other low-energy 3-body 
observables, such as the atom-dimer scattering length, can be predicted in terms of a and C. 
Near a Feshbach resonance as a — > oo, we find that C should oscillate between zero and 67.9. 
In any interval of B over which a increases by a factor of exp(7r/s ) ~ 22.7, a should have 
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FIG. 3. Contribution to a from recombination into the weakly-bound state as a function of B 
for the 907 G Feshbach resonance in 

a zero. Assuming that it is dominated by recombination into the weakly-bound state, the 
3-body loss rate should have a minimum at that value of B. If a Bose-Einstein condensate 
was prepared at such a value of the magnetic field, one could study its behavior with large 
scattering length and relatively small 3-body losses. 

We thank C.H. Greene, G.P. Lepage, J.H. Macek, and U. van Kolck for useful discussions. 
This research was supported in part by NSF grant PHY-9800964 and by DOE grants DE- 
FG02-91-ER4069 and DOE-ER-40561. 



7 



REFERENCES 



[1] Yu. Kagan, LA. Vartan'yants, and G.V. Shlyapnikov, Sov. Phys. JETP 54, 590 (1981); 

L.P.H. de Goey et al, Phys. Rev. B 34, 6183 (1986), ibid. 38, 646 (1988); H.T.C. 

Stoof et al, Phys. Rev. B 38, 11221 (1988); A.J. Moerdijk, H.M.J.M. Boesten, and B.J. 

Verhaar, Phys. Rev. A 53, 916 (1996). 
[2] P.O. Fedichev, M.W. Reynolds, and G.V. Shlyapnikov, Phys. Rev. Lett. 77, 2921 (1996). 
[3] E. Nielsen and J.H. Macek, Phys. Rev. Lett. 83, 1566 (1999). 
[4] B.D. Esry, C.H. Greene, and J.P Burke, Phys. Rev. Lett. 83, 1751 (1999). 
[5] F.A. van Abeelen and B.J. Verhaar, Phys. Rev. Lett. 83, 1550 (1999); V.A. Yurowsky 

et al, Phys. Rev. A 60, R765 (1999); E. Timmermans et al, [ |cond-mat/9805323| ] . 
[6] E. Tiesinga, B.J. Verhaar, and H.T.C. Stoof, Phys. Rev. A 47, 4114 (1993); E. Tiesinga 

et al, Phys. Rev. A 46, R1167 (1992); A.J. Moerdijk, B.J. Verhaar, and A. Axelsson, 

Phys. Rev. A 51, 4852 (1995). 
[7] S. Inouye et al, Nature(London) 392, 151 (1998); P. Courteille et al, Phys. Rev. Lett. 

81, 69 (1998); J.L. Roberts et al, Phys. Rev. Lett. 81, 5109 (1998). 
[8] J. Stenger et al, Phys. Rev. Lett. 82, 2422 (1999). 

[9] G.P. Lepage, in From Actions to Answers, TASI'89, ed. T. DeGrand and D. Toussaint, 
(World Scientific, Singapore, 1990); D.B. Kaplan, Qnucl-th/9506035| ] . 

[10] PF. Bedaque, H.-W. Hammer, and U. van Kolck, Nucl. Phys. A 646, 444 (1999); Phys. 
Rev. Lett. 82, 463 (1999). 

[11] P.F. Bedaque, H.-W. Hammer, and U. van Kolck, Nucl. Phys. A (in print) 
Dnucl-th/9906032| ] ; Phys. Rev. C 58, R641 (1998); F. Gabbiani, P.F. Bedaque, and 
HW. Griesshammer, Nucl. Phys. A (in print) [ Jnucl-th/9911034p . 

[12] K.G. Wilson, Rev. Mod. Phys. 55, 583 (1983). 



8 



